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Abstract 

In their comment to the paper "Theory of the evaporation/condensation 
transition of equilibrium droplets in finite volumes" [Physica A319, 99 
(2003)], Biskup et al. claim that in finite systems at fixed density "the 
physical significance of the conjugate thermodynamic variable is of dubi- 
ous value" . This claim is critically discussed. 
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While in the thermodynamic limit the various ensembles of statistical 
thermodynamics are trivially related to each other by Legendre transfor- 
mations, in finite systems nontrivial differences of the physical properties 
in different ensembles appear. This fact, emphasized in the preceding 
comment PQ, of course, is very well known and it is emphasized in the 
literature on computer simulations in particular, since computer 

simulations always have to deal with systems of finite size only. 

For instance, considering an Ising ferromagnet for temperatures below 
the critical temperature T c , as a function of magnetic field H, the conju- 
gate variable, the magnetization m, is a quantity that exhibits nontrivial 
fluctuations, but its average (m) increases monotonically from negative to 
positive values as H is increased from negative to positive values. Thus 
(m) = —df(T,H,L)/dH is a smooth regular function for all finite linear 
dimensions L, in the d-dimensional L d geometry with periodic boundary 
conditions, and its derivative x — {d{m) /8H)t is everywhere nonnega- 
tive and finite. When L — > oo, M(H) develops the jump singularity at 
H = from — m coex to +m coex , m coex being the value of the spontaneous 
magnetization. The details of this behavior have been elucidated in the 
literature and are well understood 0EHE1- 

The behavior of the finite Ising magnet where the magnetization m is 
held fixed as an independent variable is very different Then the con- 
jugate intensive variable H is a quantity that exhibits nontrivial fluctua- 
tions, and its average (H) = dg(T, m, L) /dm is a smooth regular function 



for all finite L, but for T <T C exhibits a loop as function of m, rather than 
increasing monotonically with increasing m. The singular behavior that 
occurs in the thermodynamic limit due to the presence of this loop was 
discussed in [SJ. Of course, if one keeps m fixed and takes L — > oo, then 
the isotherm develops a constant horizontal part H = from —m coex to 
+m C oex (This is trivially related to the corresponding result in the (H, T) 
ensemble via the Legendre transformation, of course). The different sin- 
gular behavior of H vs. m discussed in only follows in a special limit 
where one takes the limits L — ► oo and m — > — m coex (from the side where 
m > — m coex ) together, such that 5rh = (m + m CO£X )L d ^ d+V) = const. 
Since at finite L for m near ~m coex some (rounded) remnants of the 
transition present in this limit are still observed in simulations (the first 
observation was actually made in |10p. this transition and its various 
signatures are of interest for simulations, motivating the discussion pre- 
sented in In fact, in the equivalent lattice gas interpretation of the 
Ising model this transition gets the meaning of an evaporation transition 
of a droplet in a finite volume when it reaches a certain minimum size. 
The existence of this transition, first suggested in (Eqs. 26, 27), of 
that reference) was recently established rigorously |12|. 

Of course, there is nothing wrong to consider the intensive variable 
conjugate to a fixed density (magnetization density or particle number 
density, etc.) in a finite system. Considering the fluctuations of the 
temperature in a finite system in the microcanonical NVE-ensemble is 
standard textbook wisdom |13|. Similarly, the chemical potential fi in 
the NVT ensemble for finite volumes has always been of great interest 
in the context of computer simulations [3J |3J |S| and the Widom particle 
insertion method |14| has been especially devised for the sampling of the 
chemical potential in the NVT ensemble. If the statement of Biskup et 
al. that "in such systems the meaning of a conjugate variable is rather 
murky" were correct a large body of well-established work would become 
obsolete. Even if the chemical potential in a finite system has an "at best 
secondary meaning" [I] in the context of the rigorous derivations, it is a 
well accessible and useful quantity for the simulations, and it gives useful 
information for the phenomenon at hand. As an example, Fig. presents 
recent data |15| on the distribution of the chemical potential in a three- 
dimensional Lennard- Jones fluid for T ~ 0.68T C in a box of size L — 22.5a 
(a being the range of the Lennard- Jones potential) for several values of the 
particle number N near the droplet evaporation/condensation transition. 
One can see that the distribution of a changes from a single peak distri- 
bution for N = 355 representing the strongly super-saturated gas without 
a droplet to a two peak distribution near N ~ 370, where part of the time 
a droplet is present in the system, and part of the time of the sampling it 
is not present, while for N = 380 only the second peak remains, describ- 
ing the state where the droplet always coexists with the surrounding, less 
strongly supersaturated, gas. This is exactly the behavior suggested on 
the basis of the considerations described in 

In conclusion, we do not agree with the basic claim of the preceding 
comment, namely that in finite systems at fixed density "the physical 
significance of the conjugate thermodynamic variable is of dubious value" , 
but we also would like to emphasize that in science the question whether 
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something is "valuable" or not is basically subjective, and the real question 
about scientific results is whether they are right or wrong. Thus, we invite 
the reader to form his own opinion about this subject. 

Acknowledgements: I am grateful to L. G. MacDowell, P. Virnau and 
M. Miiller for a fruitful collaboration which resulted in Ref. 15 ... 
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Figure 1: Distribution p(p) of the chemical potential /i for a Lennard- Jones 
fluid at T = 0.68e/fcs and several choices of N as indicated, for a cubic box of 
volume (22.5<r) 3 , a being the range parameter of the Lennard- Jones potential 
and e its strength. From MacDowell et al. |15j . 
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